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Dijet anomaly reported by the CDF (Collider Detector at Fermilab) collaboration in 1.96 TeV pp 
collisions is investigated within the standard model by considering effects of parton distribution func- 
tions on various processes: W^+dijet, ^+dijet, WW, ZW , and top production. Since the anomalous 
peak exists in the dijet-mass region of 140 GeV with the pp center-of-mass energy y^=1.96 TeV, a 
relevant momentum fraction x of partons is roughly 0.1. In this x region, recent HERMES semi- 
inclusive charged-lepton scattering experiment indicated that the strange-quark distribution could 
be very different from a conventional one, which has been used for many years, based on opposite- 
sign dimuon measurements in neutrino-induced deep inelastic scattering. We investigated effects of 
such variations in the strange-quark distribution s{x) on the anomaly. We found that distributions 
of W-|-dijets and other process are affected by the strange-quark modifications in wide dijet-mass 
regions including the 140 GeV one. Since the CDF anomaly was observed in the shoulder region 
of the dijet-mass distribution, a slight modification of the distribution shape could explain at least 
partially the CDF excess. Therefore, it is important to consider such effects within the standard 
model for judging whether the CDF anomaly indicates new physics beyond the standard model. We 
also show modification effects of the strange-quark distribution in the LHC (Large Hadron Collider) 
kinematics, where cross sections are sensitive to a smaller- a; region of s{x). 



PACS numbers: 13.85.-t, 13.87.-a, 12.38.-t 

I. INTRODUCTION 

A signature beyond the standard model has been 
searched by the Tevatron at Fermilab (Fermi National 
Accelerator Laboratory) and by the LHC (Large Hadron 
Colhder) at CERN (European Organization for Nuclear 
Research). So far, no explicit signature has been found 
yet. However, the CDF (Collider Detector at Fermilab) 
collaboration reported a possible discovery of new phe- 
nomenon in high-pT-lepton plus dijet production with 
missing Et in pp coUisions at the center-of-mass energy 
of 1.96 TeV [l|. They showed the VF-|-dijet cross section 
and other possible contributions as a function of the di- 
jet mass nijj . Their data indicated an anomalous peak in 
the region nijj — 140 GeV. It is anomalous in the sense 
that such an excess is difficult to be explained within the 
standard model. On the other hand, such a peak cannot 
be found in another Tevatron experiment DO so that 
it may take time to settle down this issue experimentally. 

We may recollect that there was an issue of jet-cross 
section excess in a CDF experiment of 1996 Jet- 
production cross sections were measured in the pp colli- 
sions, and the data were shown as a function of jet trans- 
verse energy {Et)- Since the cross sections were signif- 
icantly larger than the QCD predications in the energy 
region Et ^ 400 GeV, the excess used to be considered 
as a signature beyond the standard model. However, it 
turned out that the jet excess could be explained if the 
gluon distribution is modified at large x (~ 0.6) where 
X is the momentum fraction of a gluon in the nucleon. 



For explaining the new CDF data of 2011, interesting 
theoretical ideas have been proposed 5] mostly by mod- 
els beyond the standard model. The ideas include new 
Z' boson, supersymmetry, Higgs boson, relation to string 
theory, and so on. However, the experience on the CDF- 
jet issue of 1996 indicates that careful studies should be 
done within the standard model, particularly by estimat- 
ing effects of parton distribution functions (PDFs). 

A fraction of the proton or antiproton beam energy 
980 GeV (=1.96 TeV /2) needs to be transferred to one 
of the dijets with the energy about 70 GeV (=140 GeV 
/2). Therefore, the relevant region of the momentum 
fraction x for initial partons is x ~ 0.1. This is a very 
rough estimate because some of the beam momentum 
should be carried, for example, by W in the final state. 
We need to check whether the PDFs are well determined 
in this X region in order to warrant the possibility that 
the CDF result is a new discovery beyond the standard 
model. 

The PDFs of the nucleon have been investigated for 
many years, and they may be considered to be established 
except for extreme kinematical conditions. In fact, the 
valence-quark distributions are rather well determined 
due to abundant charged-lepton deep-inelastic-scattering 
(DIS) data together with neutrino DIS ones by using 
strong constraints of charge and baryon-number conser- 
vations. There is a recent issue of a nuclear modification 
difference between the charged-lepton and neutrino re- 
actions [6*] . It affects the determination of valence-quark 
distributions because neutrino-nucleus DIS data are used 



for extracting "nucleonic" PDFs. However, it should not 
be large enough to affect the CDF anomaly. 

On the other hand, antiquark distributions are not well 
determined, particularly for strange and heavy quarks. 
The u and d distributions also have some uncertainties; 
however, they are relatively well determined in the region 
0.05 < X < 0.2 due to Drell-Yan experiments at Fermilab 
in combination with DIS measurements 0- Furthermore, 
a new Drell-Yan experiment is currently in progress at 
Fermilab and it could be continued at other hadron 
facilities such as the J-PARC (Japan Proton Accelera- 
tor Research Complex) 9] . The gluon distribution func- 
tion g{x) played an important role in explaining the CDF 
jet result in 1996 However, it may not be a signifi- 
cant contribution in the present case because g{x) is con- 
strained by scaling-violation measurements of the HERA 
(Hadron-Electron Ring Accelerator) except for the dis- 
tribution in the large-x region. Therefore, the first par- 
ton distribution, which we need to consider in connection 
with the new CDF anomaly, is the strange-quark distri- 
bution. 

The strange-quark distribution s{x) may have been 
thought to be determined from the opposite-sign dimuon 
measurements in neutrino- induced DIS. Current situa- 
tion is summarized in Ref. [lo|. There are some varia- 
tions among experimental groups; however, J dxx{s + s) 
is 30—50% of the light-antiquark distribution J dxx{u + 
d). By considering the neutrino measurements, all the 
global PDF analysis had been done by assuming that 
s{x) -\- s{x) has the same x-dependence with u{x) + d{x) 
until recently. There are some recent attempts to assign 
more flexible functional forms which are independent of 
u{x)-\-d{x) However, the strange-quark distribution 
is not accurately constrained as typically shown by larg e 
uncertainties of s{x) by the CTIO parametrization 

In 2008, the HERMES collaboration cast doubt on the 
conventional strange-quark distribution [l2| . They mea- 
sured semi-inclusive hadron production in charged-lepton 
DIS and they obtained a very different strange-quark 
distribution from the ones by global analyses with the 
neutrino DIS data. The distribution s{x) is much softer 
than usual global- analysis distributions at (5^=2.5 GeV^. 
For example, the s{x) is two-times larger than the usual 
one at x = 0.03 and it is much smaller at x = 0.2. Of 
course, hadron-production analysis depend much on the 
employed fragmentation functions, which have large un- 
certainties for kaon functions [isj . Such an ambiguity 
needs to be considered in evaluating the HERMES result 
seriously. 

The light antiquark distributions u{x)^ d{x), and s{x) 
are considered to be equal in perturbative QCD (pQCD) 
because they are produced from the gluon splitting pro- 
cess g ^ qq and their masses are expected to be negligible 
in comparison with a hard scale of a reaction. Of course, 
higher-order pQCD corrections could give rise to differ- 
ences between u, d, and s because of evolution with a 
mixing term q ^ q in a splitting function at the next-to- 
leading order (NLO) together with differences between 



the initial distributions u, d, and s Q. However, such 
NLO effects are not large enough to explain the mea- 
sured differences, so that a non-perturbative mechanism 
is needed for explaining the differences between u, d, and 
s. For the details of these non-perturbative models, the 
reader may look at Ref. 7]. 

It is the purpose of this article to discuss whether varia- 
tions of s{x), such as the one suggested by the HERMES, 
could affect the CDF anomaly. In particular, cross sec- 
tions of VF-l-dijet, Z-l-dijet, WW, ZW, and top produc- 
tion processes are calculated by using the code GR@PPA 
(GRace At Proton-Proton/ Antiproton collisions) [l3| by 
taking three functional forms of s{x). The GR@PPA is 
an event generator for V {W and Z)-|-jets, VV, and QCD 
multi-jet production processes in pp and pp collisions. 

This article is organized in the following way. In Sec. 
im our method is explained for calculating the cross sec- 
tions by the GR@PPA, and employed strange-quark dis- 
tributions are explained in Sec. IIIII Results for the 
lepton-plus-dijet cross sections are shown in Sec. llVi and 
they are summarized in Sec. |Vl 



II. METHOD FOR CALCULATING CROSS 
SECTIONS 

In describing the cross sections for M^-|-dijet, Z-|-dijet, 
WW, ZW, and top productions, the event generator 
GR@PPA [IJ] is used in this work. First, this event gen- 
erator is explained in Sec. Ill Al We use the CTEQ6L1 
distributions for the PDFs in the proton and anti-proton. 
Next, in order to find effects of the strange-quark dis- 
tribution, different forms of s{x) are provided in the 
GR@PPA code as explained in Sec. Hm The other PDFs 
are not modified from the CTEQ ones. 



A. Event generator GRQPPA 

The GR@PPA is a framework to calculate cross sec- 
tions of proton-proton/antiproton collisions by using the 
GRACE [l^, in which Feynman diagrams are automat- 
ically computed at the one-loop level as well as the tree 
level in the standard model. The GRACE is intended 
mainly for lepton collisions, so that it is implemented by 
including features of hadronic collisions such as the PDFs 
in the GR@PPA [3. 

The cross section for a certain hadron-hadron reaction 
process can be calculated by the sum of partonic subpro- 
cess contributions 




(1) 

where fi'{xh,Q'^) is the PDF of the hadron h, Xh is the 
momentum fraction carried by the parton i or j, is 
a hard scale in the reaction, and daij^Fis)/d^F is the 
cross section for the partonic subprocess i+j — >■ F. Here, 
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i and j indicate parton types, F is the final state, s is 
the square of the center-of-mass energy in the parton level 
(s — {pi + pj)'^), and is the phase space of the final 
state. 

As for the partons, we take into account If types (u, 
d, s, c, b, g, u, d, s, c, and b). There are a huge num- 
ber of subprocesses which contribute to the cross sec- 
tion. It is obvious that some of them have identical dia- 
grammatic structure with slight changes in coupling con- 
stants such as the CKM (Cabibbo-Kobayashi-Masukawa) 
factors. Therefore, computation time is significantly re- 
duced by calculating selected base-subprocesses and then 
weight factors are simply multiplied: 

^ /■ , , dalf^T^is;m,a) 
fT = > / dxi ax2 a<Pp WijF ■ (2) 

Here, daf°^j^p'^'^{s; m, as) is the cross section for the base 
subprocess at the cm. energy squared, masses, and cou- 
plings. The Wijp is the weight factor which takes into 
account the PDFs, CKM matrix Vckm, and decay fac- 
tors of W/Z bosons: 

Wi,F = fUxuQ^) ff{x2,Q^)\VcKMf'' , (3) 

where K is the number of W boson in the process. 

The GR@PPA code can be obtained from its web page 
[l^ . Outline of GR@PPA is explained in the following; 
however, rnuch details should be found in the original 
articles [13, [T3| • A cross section can be calculated in the 
following steps. 

(1) Particles and vertices in theory 

First, a theoretical model is supplied in the form of 
particle content, parameters, and interaction ver- 
tices derived from a model Lagrangian. In this 
work, the model is the standard model of strong 
and electro-weak interactions. The model informa- 
tion is supplied as model files. 

(2) Input for particles and kinematics 

Initial- and final-state particles as well as reaction 
kinematics need to be supplied for calculating cross 
sections. The kinematical conditions are cm. en- 
ergy, kinematical cuts for jets and final-state par- 
ticles, hard scales for the running coupling con- 
stant and PDFs, renormalization and factorization 
scales, and the others. They are supplied in the files 
grcpar.F and upinit.F as instructed in Table II of 
Ref. The kinematical cuts should be supplied 
as they were used in the CDF analysis. 

(3) Diagram generation 

From the supplied model information on the parti- 
cles and vertices, Feynman diagrams are automat- 
ically generated by using mathematical techniques 
of graph theory. Here, a node corresponds to a ver- 
tex or an external particle, and an edge does to a 
propagator or a connection between a vertex and 
an external particle. First, vertices are generated 



for a given condition of a process. Then, the ver- 
tices are connected by considering all the possible 
processes to satisfy obvious conservations such as 
electric charge and fermion number. Finally, parti- 
cles are assigned to the propagators. 

(4) Matrix-element generation 

The generated Feynman diagrams are automati- 
cally calculated by using spinors and gamma ma- 
trices in a numerical way. Based on the Feynman- 
diagram configuration, the code calls the model li- 
braries in (1) for calculating each amplitude. Here, 
propagators are expressed by bi-linear forms of 
wave functions, Dirac spinors for fermions, and po- 
larization vectors for spin-1 bosons. 

(5) Monte Carlo integrals 

Finally, the phase-space integrals of Eq. ([2]) 
have been calculated by the Monte-Carlo method 
BASES. The parameter NCALL is the number of 
sampling points in integration grids, and it should 
be large enough to obtain convergent numerical re- 
sults. The recommended values of NCALL are sup- 
plied in the subroutine grcpar.F, but they may be 
changed according to kinematical conditions. 

In the published GR@PPA code, the single-top pro- 
duction is not supplied, so that a separate code is created. 
In this work, the single-top and tt production processes 
are calculated, and resulting distributions are shown as 
top events. 

Since there are many processes which contribute to the 
p+p — > {ii^) + 2 jets -I- we cannot show all of them 
here. Only typical processes are shown. For example, one 
of M^-|-dijet processes is shown in Fig. [1] where W~ is 
created in an intermediate stage to produce an electron 
and an anti-neutrino in the final state, and two jets are 
produced by a qq pair from a hard gluon. Several more 
processes arc shown in Appendix \X\ for providing an idea 
of other contributions. For getting information on all the 
processes, one may look at figure files created by running 
the GR@PPA. In the same way, typical contributions of 
Z+dijet, top, WW, and ZW processes are shown in Figs. 
H nil and O respectively. 

In our work, elementary partonic cross sections are cal- 
culated by the GR@PPA. Subsequent parton branching 
and final fragmentation into hadrons are not included. 
They could be calculated, for example, by using the event 
generator PYTHIA [l^ . It is the purpose of this work to 
show gross properties of strange-quark effects simply by 
calculating the hard process part. In future, if much de- 
tailed comparisons become necessary with the CDF and 
DO experimental data, they needs to be considered. This 
work is the first step toward such an approach. 

B. Applied kinematical conditions in GR®PPA 

In obtaining the cross section, various kinematical cuts 
are applied. They are the same as the ones used in the 
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FIG. 1: Typical Vy+dijet process. The notation q indicates a 
quark u, d, s, or c. 
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FIG. 2: Typical .Z+dijet process. 




FIG. 3; Typical top-production process. 




FIG. 4: Typical WW process. 
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FIG. 5: Typical ZW process. 



CDF analysis [T| . In order to search for the signal p + 
p — ^ WW/ZW — >■ ivg + dijets, events with (two energetic 
jets) + (one high-py electron or muon) + (missing Et) 
are selected. The high-p^ lepton from W needs to be 
accompanied by missing transverse energy (fx ) for a 
non-interacting neutrino. 

The momentum cut p^ > 20 GeV is applied for re- 
moving QCD backgrounds from semi-leptonic decays of 
hadrons and jry^j < 1.0 is due to detector ranges. For se- 
lecting the leptonic W candidate, a cut is also applied for 
the missing Et as I^t > 25 GeV in order to remove the 
QCD backgrounds. Furthermore, the transverse-mass 
cut > 30 GeV also reduces the QCD backgrounds. 
The cone radius AR ^ v^(A0)2 ^ {Arj^ = 0.4 is used for 
the jet definition with the cuts for each jet as Ei^*' > 30 
GeV and rjj < 2.4. Additional cut conditions are ap- 
plied as lAijjjl < 2.5, A<l)jf^j^ > 0.4, ARjg > 0.52, 
and PTjj > 40 GeV for rejecting back-to-back events, 
multi-jet (fake missing I^t) events, semi-leptonic decays 
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of hadrons, and for better agreement between data and 
simulation, respectively. Gaussian smearing effects are 
taken into account for W and Z by assigning an effective 
width of 15 GeV in the dijet-mass distribution. 

III. STRANGE-QUARK DISTRIBUTION 
FUNCTION 

A. Definition and experimental information 

The quark distribution function qi{x) with flavor i is 
defined by the amplitude for the nucleon (N) to emit the 
quark i and then to absorb it at a separated light-cone 
coordinate y~ [uj 

X7+p^,(0+,0-,0i)|7V(p)), (4) 

where p is the path-ordered (P) link to make the expres- 
sion gauge invariant: 

p = Pexpi^-igJ^ dz-A+(0+,z-,Oi)ta|. (5) 

The indicates the light-cone coordinates — (a° ± 
a'^)/\/2, and a± is the two-dimensional transverse coor- 
dinate. The ip and A are quark and gluon fields, respec- 
tively, and ta is the generator of SU{3) group and it is 
given by the Gell-Mann matrix Aa as ta = Aa/2. 

The PDFs can be calculated in nonperturbative theo- 
retical methods such as lattice QCD, bag model, chiral 
soliton model, and so on. However, obtained PDFs are 
not accurate enough to supply precise cross sections, so 
that global analyses of world data have been done for 
obtaining the optimum PDFs. 
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FIG. 6: Strange-quark distribution in neutrino-induced 
opposite-sign dimuon production process {Vfj, -|-p — >■ -I- 
X). 



As explained in Sec. U the important PDF to be in- 
vestigated is the strange-quark distribution s{x) in con- 
nection with the anomalous CDF result. The strange 
distribution has been determined from dimuon measure- 
ments in neutrino DIS (i^^ + p ^ A^^A'^ + ^) [IS 
shown in Fig. |6l The dimuon production occurs from an 
initial strange (or down) quark to produce an interme- 
diate charm quark which decays into and s (or d). 
From measurements, the strange-quark-moment ratio to 
the light-antiquark one 

^ Jdxx[six,Q^) + six,Q^)] 

Jdxx[u{x,Q2) + d{x,Q^)Y ^ ' 

was obtained. This ratio varies depending on experi- 
ments, but it is in the range of 0.3 ^ k ^ 0.5 p^ . 

This neutrino process has been the only constraint for 
determining the strange distribution s{x) in global PDF 
analyses. However, the HERMES collaboration recently 
indicated that the distribution s{x) could be much dif- 
ferent from the current one (l2j by semi-inclusive kaon 
production in Fig. [71 The shaded blob in the figure indi- 
cates a fragmentation process, which is described by the 
fragmentation functions (FFs). It should be noted that 
the uncertainties of the FFs are very large for the kaon 
as shown in Ref. [l3l |. For example, such an uncertainty 
effect is illustrated in determining the polarized strange- 
quark distribution As (a;) from data sets including semi- 
inclusive measurements [ISJ. In any case, considering the 
large differences between s{x) of the HERMES and the 
conventional one, we think that the strange distribution 
could have much larger uncertainties than the ones sug- 
gested by most PDF analyses [Tl| . 

B. Modifications of strange-quark distribution 

We discuss modifications of the strange-quark distri- 
bution s{x) -\- s{x) obtained in a typical global analysis. 
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FIG. 7: Semi-inclusive kaon production in charged-lepton 
deep inelastic scattering {e + p ^ e' + K + X) . 
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We define the weight function or modification function 
Wi{x,Q^) for the strange-quark distribution by 



[s{x,Q^) + s{x,Q% 

= Wi{x, Q^) [s{x, Q^) + s{x, Q^)]pDF set , 



(7) 



at any Q^, where perturbative QCD can be applied. 
Here, the CTEQ6L1 distributions [l^ are used for the 
PDF set of the global analysis because they have been 
used in the GR@PPA code. Namely, the function Wi in- 
dicates a modification from the s + s distribution of the 
CTEQ global analysis. The index i indicates a modifica- 
tion type (i =1 or 2) explained in the following. 

As the first possible choice, the weight function is de- 
termined from measurements of the strange-quark distri- 
bution s{x) + s{x) by the HERMES collaboration. The 
average of the HERMES data is 2.5 GeV^, so that 
we define Qq = 2.5 GeV^ as the scale in providing the 
initial PDFs including the modification. The follow- 
ing hyperbolic-tangent function wi is assigned for the 
modification by looking at the differences between the 
CTEQ6L1 distribution and the HERMES data in Fig. E 



wi {x, Qq) = 1 — tanh 



Xq 



(8) 



The HERMES a;(s -I- s) distribution is much larger than 
the CTEQ6L1 at x < 0.06 and it is much smaller at 
X > 0.09. The parameter xq indicates the cross-over 
point, and Ax does the width of the modification curve. 
These two parameters are determined so as to fit the 
HERMES data, and obtained values are listed in Table 
H The obtained weight function wi{x,Qq) is shown in 
Fig. in and it explain the data well except for slight 
deviations at x ^ 0.3. 

In addition to the function wi, we consider the follow- 
ing one: 



W2{x,Ql) 



xo 



1 -I- tanh , ^ 
Ax 

[u{x, Ql) + d{x, Qo)]PDF sot 



[s{x,Q 



PDF set 



(9) 



This function ^2(2;) has the opposite functional form to 
wi{x) in the sense that the small-x part is suppressed and 
the large-x distribution is enhanced to reach to the u + d 
distribution. This functions W2 ■ x{s + s) is shown in Fig. 
[5] in comparison with the HERMES data, wi ■ x{s + s), 
and x{u + d). Although W2{x,Qq) is not supported by 



TABLE I: Parameter values in the weight functions of Eqs. 
dH), (O, and Uni) at Q^=2.5 GeV^ and = 2M^. 



Function 


xo 


Ax 


CO 


Cl 


wi{x,Qi) 
W2(a;,Qo) 

W2{x,Q^) 


0.0796 
0.0796 
0.0768 
0.0727 


0.0253 
0.0253 
0.0408 
0.0220 


0.802 
1.589 


-0.475 
0.417 




FIG. 8: x{s + s) distribution functions are compared with 
HERMES dataM and x{u + d). The CTEQ6L1 distribu- 
tion x{s + s) [19 1 is shown by the solid curve, and the mod- 
ified distributions by wi{x,Qo) and W2{x,Qq) are shown by 
the dashed and dash-dotted curves, respectively. The weight 
functions wi{x,Qa) and ■W2{x,Qq) are explained in the main 
text. For comparison, the x{u + d) distribution is shown by 
the dotted curve. 



the HERMES data, we simply consider it as a trial func- 
tion to investigate sensitivity of the CDF result on the 
strange distribution. The strange distribution is roughly 
constrained by the neutrino measurements; however, its 
x-dependent functional form is not determined. In par- 
ticular, the large- a; region is not explored experimentally. 
As there is an issue of a significant intrinsic charm distri- 
bution at large x 0, [2lj , there may be an enhancement 
in the strange at large x From Fig. [U the small- a; 
distribution may seem to be extraordinarily small. How- 
ever, if the function W2 is evolved to the average (^20 
GeV^) of CCFR and NuTeV measurements, the suppres- 
sion is not very large as it becomes clear in Sec. IIII CI 
It should be also noted that these functions wi and W2 
roughly correspond to the bounds of the s{x) uncertainty 
in the CTIO distribution [11,] . 



C. Q evolution 

The hard scale in providing the PDFs for the CDF 
measurements is much larger than Qq — 2.5 GeV^, where 
the HERMES data are provided for the s + s distribu- 
tion, so that the strange-quark distribution needs to be 
evolved. As a typical hard scale in calculating the CDF 
dijet cross sections, we take Q^ — 2M^ by considering 
the dijet mass rrijj ^ 100 GeV and the boson (W and 
Z) masses. The standard DGLAP (Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi) equations are used for the Q^ 
evolution by using the code of Ref. [2^. In the case 
without the modification {w — 1), we checked that the 
evolved distribution is numerically consistent with the 
CTEQ one. The evolved strange distributions are shown 
in Fig. ini The strange quarks are copiously produced 
especially at small x due to Q^ evolution, so that the 
distributions become large at small x. Therefore, the 
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FIG. 9: evolved x{s + s) distribution functions are shown. 
The CTEQ6L1 distribution [l^ and the modified distribu- 
tions by w\{x,Q'q) and W2{x,Q'^) are evolved to = 2M^. 
For comparison, the xiu -\- d) distribution is shown at — 
2M^ by the dotted curve. 



ments although one should note that the neutrino- 
experimental ratio of Eq. ([6]) is for the second moments. 
At = 2M^, the modifications become much smaller 
and they are in the range 0.8 ^ wi^2 ^ 1-2 at small x 
(< 0.08). The modifications also become smaller with 
increasing at large x (> 0.1); however, variations 
are not as large as the ones at small x. 

For including the modifications into the event genera- 
tor, the evolved weight functions are fitted by the func- 
tions 



Wi{x,i 
W2{x,i 



-cq — ci tanh 



Co 



ci tanh 



X — Xq 

Ax 

X - 



Xq 



Ax 

[u{x, Q^) + d{x, Q^)]pDF set 

[s{x, Q^) + s{x, Q^)]pDF set '' 



(10) 



differences between (s + s)pdf set, (s + s)i, and (s -t- 5)2 
become smaller than the original ones at Qq . In compar- 
ison, the large-x (> 0.1) distributions are not so signifi- 
cantly modified by the evolution. 

In order to clearly show the evolution of the weight 
functions, we divide the evolved strange distributions 
Si{x,Q'^) by the CTEQLl distribution to calculate the 
weight function Wi{x,Q'^) as given in Eq. ([7]) at any 
Q^. In Fig. [TUl three values are taken for show- 
ing the evolution of the weight functions wi and 
W2. They are Q^=2.5 GeV^, 20 GeV^, and 2M^ for 
the initial HERMES scale, NuTeV scale, and CDF-jet 
scale, respectively. Although the modifications are large 
at (5^=2.5 GeV^ and it is in the range < wi,2 ^ 2, 
they become smaller as increases. At the neutrino 
DIS scale of = 20 GeV^, the modifications becomes 
smaller (0.5 < W1.2 1-5) at small x (< 0.08). In com- 
parison with u -I- d at = 20 GeV^, the strange to 
light-antiquark ratio is 0.23 ^ {s + s)/{u + d) < 0.68, 
which is not very different from the neutrino measure- 
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FIG. 10: The weight functions wi{x,Q^) and W2{x,Q'^) are 
shown at Q^=2.5 GeV^ 20 GeV^ and 2M^ by the solid, 
dotted, and dashed curves, respectively. For comparison, the 
ratios (s + s)/{u + d) are shown at Q^=2.5 GeV^ and 2M^. 



where cq, ci, Xq, and Ax are the parameters determined 
by the fitting. Obtained parameters are listed in Table 
m Since it takes time to calculate the cross sections by 
the Monte Carlo integrals in the GR@PPA, these fitted 
functions are used instead of calculating the evolution 
within the Monte Carlo code. 



D. Comments on recent PDFs 
and theoretical models 

We comment on the strange-quark distribution in re- 
cent global analyses. Until a few years ago, the dis- 
tribution s{x) is simply assumed to be proportional to 
[u{x)+d{x)\/2 together with the assumption s{x) = s{x). 
However, more flexible parametrizations were recently in- 
vestigated in some analyses by allowing independent x- 
dependence for s{x). Then, the uncertainty of the de- 
termined s{x) became large in general. The size of the 
uncertainty of s{x) depends on analysis groups. A large 
uncertainty was suggested in the CTIO PDF set for the 
strange-quark distribution as shown in Fig. 12 of Ref. 
[ll| although uncertainties of other analysis groups are 
smaller. The CTIO indicates a large uncertainty range 
< w{x, Q^) < 2 at the scale = 4 GeV^ depending 
on the X region. This uncertainty range roughly corre- 
sponds to the functions wi^2{x, Q^) in Fig. [101 Therefore, 
although the modifications seem to be very large espe- 
cially at =2.5 GeV^ in Figs. [8land[T0l they are allowed 
by the current experimental data. As it is explained in 
recent global- analysis articles [llj, that Z and W produc- 
tion cross sections at LHC are sensitive to the strange- 
quark distribution, so that this issue may be solved in 
future LHC measurements [23 |. 

Theoretically, the strange-quark distribution stems 
from perturbative and non-perturbative mechanisms. 
The perturbative creation of ss pairs is described by the 
DGLAP evolution equations. The non-perturbative 
contributions are calculated by phenomenological models 
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such as meson clouds, chiral soliton, and so on; how- 
ever, obtained distributions inevitably depend on the 
used model. The lattice QCD method has been mak- 
ing progress; however, it is not at the stage to provide 
the detailed x dependent functional form. A few low- 
moments of the PDFs have been obtained so far in the 
lattice QCD [25^]. 

As a nonperturbative effect on the PDFs, an intrinsic 
charm distribution has been often discussed. The charm 
distribution c{x) is usually considered to be created by 
a perturbative mechanism by the gluon splitting g 
cc; however, there could exist the intrinsic charm in the 
proton according to the light-cone Fock space picture [l^l 

I p >= I uud > + • • • +1 uudcc >+■■■. (11) 

In addition, a meson-cloud picture also suggests the in- 
trinsic charm distribution by p(uud) — ^ D^(uc)K^{udc), 
p{uud)J/ip{cc) [l^. In the same way, the nonperturba- 
tive strange-quark distribution could exist by considering 
the strange-quark content in the proton as [24I 

I p >= I uud > + • • • +1 uudss >-!-•••. (12) 

It is known that the intrinsic charm distribution appears 
as a bump in the large-x region as shown in Figs. 2 and 
3 of Ref. [2l|. In the same way, the strange-quark dis- 
tribution could have a distribution in the large-x region, 
which may be related to the CDF anomaly. These non- 
perturbative studies have been done for light antiquark 
{u and d) distributions 0, [l^] and for the strange distri- 
bution [2^ as well as for the charm In this article, 
we do not step into the details of these mechanisms and 
their comparisons with the CDF anomaly. 

In future, it is important to investigate the details 
of the theoretical models and the lattice QCD for the 
strange-quark distribution, and then their connection to 
the CDF anomaly. Since there is no established theory 
for the strange-quark distribution, we simply vary the 
CTEQ strange distribution to the softer one of the HER- 
MES collaboration and to the harder one as explained in 
Sec. ImBl 

IV. LEPTON-PLUS-DIJET CROSS SECTIONS 

Using the GRQPPA code and the modified strange- 
quark distributions, we calculate each process contribu- 
tion to the lepton+dijet cross section. In particular, 
W-l-dijet, Z-l-dijet, top, WW, and ZW cross sections 
are calculated. Using the CTEQ6L1 PDFs without the 
strange modification, we obtain the events as a function 
of the dijet mass in Fig. [TT] We checked that numerical 
results are consistent with other ones, for example, by us- 
ing the ALP GEN [28*1 by calculating these cross sections 
for Tevatron and LHC kinematics. 

The lepton-plus-dijet cross sections are dominated by 
the Ty-f dijet processes and other contributions are small. 
The ordinate indicates the event rate in the dijet mass 



1500 




(GeV) 

FIG. 11: Contributions to lepton-plus-dijet events in the CDF 
kinematics. The W+dijet, Z+dijet, top, WW, and ZW 
events and their total are shown. Here, the default strange- 
quark distribution is used without modification. 



interval of 8 GeV. Since the detector acceptance informa- 
tion is not available for public, and since parton shower 
and final fragmentations are not included in our calcula- 
tions, the magnitude of the ordinate cannot be compared 
with the CDF measurements. Therefore, we may look at 
overall shapes and each contribution. We find that the 
overall features agree with the CDF measurements. 

Next, we show the strange-quark effects on the cross 
section. The evolved weight functions wi{x,Q'^) and 
W2{x, Q^) are used for calculating the strange-quark dis- 
tributions in the GR@PPA. Obtained event rates are 
shown in Fig. [T^l It indicates that the cross section 
increases if the hard strange-quark distribution S2{x, Q^) 
is used, and it decreases if the soft distribution si(x, Q^) 
is used. As mentioned in the introduction section, the x 
region which affects the CDF result is x ^ 0.1. As shown 
in Fig. ini the distribution s{x) significantly increases in 
the hard strange distribution S2{x) and decreases in the 
soft distribution si{x) at a; ^ 0.1. This fact gives rise to 
the modifications in the lepton+dijet cross sections. 

ISOOt 
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mj.(GeV) 

FIG. 12: Effects of the strange-quark distribution on lepton- 
plus-dijet events in the CDF kinematics. Three curves are 
shown by taking the strange-quark distribution s{x) without 
modification and the distributions multiplied by ■wi{x,Q'^) 
and W2{x, Q^). 
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FIG. 13: Effects of the strange-quark distribution on VF+dijet 
events in the CDF kinematics. Three curves are shown by tak- 
ing the strange-quark distribution s{x) without modification 
and the distributions multiplied by wi{x,Q^) and W2{x,Q'^). 
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FIG. 14: Effects of the strange-quark distribution on W^-|-dijet 
events in the LHC kinematics. The notations are the same in 
Fig. [131 



It is interesting to find that the distribution shape as 
the function of the dijet mass becomes wider if the hard 
strange-quark distribution is used. It has a tendency to 
partially explain the anomalous CDF excess, although a 
sharp CDF-like bump is difficult to be obtained within 
modifications of PDFs. However, because the CDF find- 
ing is in the shoulder region of the cross section, a slight 
shift in the cross-section shape may explain the CDF ex- 
cess if they accumulate enough data to obtain accurate 
cross sections. 

We mentioned a possible bump structure of s{x) in 
Sec. nil Dl According to our experience on the numerical 
estimates, it is not easy to explain the sharp CDF peak 
by a bump in s{x) at medium x, even if it exists, because 
cross sections are obtained by integrating the PDFs over 
a certain region of x, more specifically over a certain ra- 
pidity range. In addition, the evolution smears out 
the bump. 

The modification effects are calculated for each pro- 
cess. We found that there are sizable effects on the 
W+dijet and Z-|-dijet cross sections from the strange- 
quark modification. However, the effects are very small 
in top, WW, and ZW production processes. Since the 
magnitude of the cross section is largest in W^-|-dijet, we 
show the modification effect on this process in Fig. [T^l 
Because the lepton-|-dijet cross section is dominated by 
the VF-l-dijet process, the overall modification effects are 
similar to the ones in the total event rates in Fig. [T^ 
The hard strange-quark distribution increases the cross 
section and shifts the shape toward the larger dijet-mass 
region, whereas the soft distribution has opposite effects. 

Next, the same W^-l-dijet cross sections are calculated 
for the LHC kinematics of 14 TeV cm. energy. Two 
replacements in the cm. energy and initial state, ^/s — 
1.96 TeV — )■ 14 TeV and pp — >■ pp, are simply done in the 
GR@PPA code for the calculation. The other cut con- 
ditions are kept the same simply for finding kinematical 
effects. The obtained results are shown in Fig. [TH We 
notice large differences from the the event rates of the 
CDF kinematics in Fig. [T31 The cross section decreases 



(increases) by the hard (soft) strange-quark distribution 
and the shape becomes narrower (wider). These tenden- 
cies are opposite to the CDF ones in Fig. [T31 This is sim- 
ply caused by kinematical effects. Because the cm. en- 
ergy is much larger in the LHC, the process is sensitive to 
the smaller-x region, particularly at x ^ 0.02. Here, the 
hard distribution S2{x) becomes smaller than the stan- 
dard one without the modification at a; ~ 0.02, and the 
soft one si{x) becomes larger as shown in Fig. [SI These 
differences produce the opposite strange-modification ef- 
fects to the ones of the CDF kinematics. 

These results indicates an interesting conclusion that 
the CDF analysis should be carefully done by considering 
the uncertainty of the strange-quark distribution. Fur- 
thermore, the strange-quark distribution can be tested 
by future LHC experiments. In particular, the extraor- 
dinary HERMES result on s{x) can be tested by both 
Tevatron and LHC measurements on vector-boson pro- 
ductions [2J|. It is important to note that Tevatron and 
LHC probe different x regions of the strange-quark distri- 
bution. There is still a possibility that the CDF anomaly 
could be explained within the standard model in terms 
of the PDF modification if the peak in the dijet-mass 
spectrum is not too sharp. 

We did not step into the actual calculations based on 
some theoretical models. We simply explored effects due 
to the modifications of the strange-quark distribution on 
the lepton-|-dijet cross sections. In future, it should be 
interesting to theoretically investigate the strange-quark 
momentum distribution by hadron models and possibly 
by lattice QCD. We leave them for our future projects. 



V. SUMMARY 

We investigated a possibility that the anomalous CDF 
result in the lepton-t-dijet cross section can be explained 
within the standard model by modifications of the par- 
ton distribution functions. In particular, the effects of the 
strange-quark distribution have been calculated by using 
the event generator GR@PPA. The strange-quark distri- 
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bution s{x) has been thought to be determined by the 
opposite-sign dimuon measurements in neutrino deep- 
inelastic scattering. However, it became obvious by the 
HERMES collaboration that its momentum distribution, 
namely the x dependence, has not been reliably deter- 
mined at this stage because the much softer distribution 
was suggested by their data. 

In this work, two modifications were considered for the 
strange-quark distribution. One is the soft distribution 
suggested by the HERMES collaboration, and another is 
the hard distribution. The modified distributions were 
evolved to the ones at the scale = 2M^. The hard 
distribution increases the lepton-|-dijet cross section and 
also spread the cross-section shape toward the large dijet- 
mass region. It tends to partially explain the CDF ex- 
cess. The sizable strange-quark-modification effects were 
obtained for the subprocesses of VF-|-dijet and Z+dijet 
productions, whereas the effects are very small for the 
subprocesses of top, WW, and ZW productions. We 
found that the PDFs in the region of x ~ 0.1 affect the 
CDF lepton-t-dijet cross sections. 

The modification effects were also calculated for the 
LHC kinematics, and the interesting results were ob- 
tained in the sense that opposite effects are found to the 
Tevatron case. The hard distribution decreases the cross 
section and shrinks the dijet-mass spectrum, which is op- 
posite to the ones in the Tevatron kinematics. This is 
caused by the fact that the PDFs in the region x ~ 0.02 
affect the lepton-f dijet cross section in the LHC kinemat- 
ics. These studies indicate that the x dependence of the 
strange-quark distribution can be investigated by the W 
and Z production measurements in both Tevatron and 
LHC. These two- facility experiments probe different x 
regions of the strange-quark distribution. 

There are still large uncertainties originating from the 
strange-quark distribution. If the CDF anomalous peak 
is not too sharp, there is still a room to explain the excess 
within the standard model in terms of the PDFs without 
introducing any new physics beyond the standard model. 
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Appendix A: An example of Vy+dijet processes 

There are too many processes, which contribute to the 
W-|-dijet cross section, to show all of them. The figure 
files to show all the processes are automatically generated 
by running the GR@PPA code. The reader may look at 
the processes in postscript files created by running the 
GR@PPA. One of the files is shown in Fig. [HI 




FIG. 15: Some of W-|-dijet processes are shown by taking a 
figure file created by running the GR@PPA code. These are 
just a part of all the calculated processes for the M^-|-dijet 
events. 
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